function f= curvefun( x,t )
    f=x(1)+x(2)*exp(-0.02*x(3)*t);
end

